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Iron-chromium forms the basis of most of the stainless steel grades in the markets. Recently new insights 
into the physical and chemical properties of Fe-Cr based alloys have been obtained. Some of the new results are 
quite unexpected and call for further investigations. The present study addresses the magnetic contribution in 
the atomic driving forces behind the chemical composition in Fe-Cr alloyed with Al, Ti, V, Mn, Co, Ni, and Mo. 
Using the ah initio exact muffin-tin orbitals method and an Ising-type spin model, it is found that the magnetic 
moment of the solute atom combined with the induced changes in the magnetic moments of the host atoms form 
the main framework in determining the mixing energy and chemical potentials of low-Cr Fe-Cr based alloys. 
The results obtained in the present work are related to tuning of the microstructure and corrosion protection of 
low-Cr steels. 



I. INTRODUCTION 

The superb properties of stainless steel, like resistance to 
corrosion, high strength, ductility, low maintenance, and rel- 
atively low cost make it an ideal base material for a host of 
commercial applications. For example, stainless steels are 
used in cook-ware, cutlery, hardware, surgical instruments, 
major appliances, industrial equipment, as an automotive and 
aerospace structural alloy, and construction material in large 
buildings. 

Stainless steel differs from corrodible carbon steel mainly 
by the amount of chromium in the bulk. Stainless steels have 
sufficient amount of chromium present so that a thin and trans- 
parent film of chromium oxide rapidly forms to the open sur- 
face of the metal which prevents further surface corrosion and 
blocks corrosion from spreading into the internal structure of 
the material. More importantly, this oxide layer quickly re- 
forms when the surface is scratched. This phenomenon is 
called 'self-healing passivation' and is seen also in other met- 
als, such as aluminum and titanium. 

High oxidation-resistance in air at ambient temperature is 
normally achieved with chromium additions in steels. The 



observed steep increase of the corrosion resistance of the fer- 
ritic stainless steels starts when the Cr content in bulk reaches 
the level of about 10 al.%A This bulk threshold of Cr corre- 
lates surprisingly well with the calculated reversal point of the 
magnitudes of the bulk and surface effective chemical poten- 
tial differences. This reversal of the chemical potentials is the 
cause for the outburst of Cr on the otherwise pure Fe surface of 
Fe-Cr alloys^ In addition to Cr many other elements are fre- 
quently used as minor alloying additions in stainless steels to 
further improve the properties of stainless steels grades. How- 
ever their impact on the chemical potential of the host matrix 
is not well documented. 

Alloying elements in steels can be divided into two main 
categories namely austenite and ferrite stabilizers. From the 
elements discussed in the present work Cr, Ti, Mo, V, and Al 
belong to the ferrite stabilizers and Ni, Mn, and Co are austen- 
ite stabilizers. In addition to the effect on the austenite/ferrite 
stabilization, these elements affect many other properties of 
steels. In the following we list some typical applications of al- 
loying. Aluminum is often used in high temperature corrosion 
resistant materials 3-5 due to its ability to form a highly stable 
and protective oxide scale on the open surface when exposed 



2 



to oxidizing environment. The physical properties and corro- 
sion resistance of Fe-Cr-Al as a function of the chemical com- 
position of the alloy have been studied quite extensively^— 
Titanium is an element added to steels because it increases the 
strength and resistance to corrosion. Furthermore, Ti provides 
a desirable property to alloys: lightness. Its density is less 
than half that of steel, so a titanium-steel alloy weighs less 
than pure steel and is more durable and stronger. These physi- 
cal and mechanical properties make titanium steels a possible 
material for example in aircraft and spacecraft engineering, 
chemical industry, and medicine* 2 ^— Vanadium is a com- 
mon component, for instance, in tool steels. Vanadium ad- 
dition to Fe-Cr increases the hardness, tensile strength, wear 
resistance, and impact resistancei 2 ^ 2 ^ However, little has been 
published on the corrosion properties of the Fe-Cr-V alloys 
and on the role of vanadium in the passivation process^— 
Molybdenum is used to increase the strength to weight ratio 
and weldability^ 

Nickel is the most frequent alloying element to stabilize 
Fe-Cr alloys in the austenite structure. This crystal structure 
makes such steels paramagnetic and less brittle at ambient 
conditions. Also manganese has been used in many stain- 
less steel compositions. Manganese preserves the austenite 
structure in steels as does nickel, but at a lower cost. Austen- 
ite grades are perhaps the most widely used stainless steels in 
modern applications. 23,34 Cobalt is used to increase heat and 
wear resistance (improve anti-galling property) and to pro- 
duce excellent magnetic properties and good ductility— 

In this Report, we investigate the role of minor alloying el- 
ements (Al, Ti, V, Mn, Co, Ni, Mo) on the relative energetics 
of Fe and Cr and its possible consequences to surface prop- 
erties. The basic data is obtained by an ab initio electronic 
structure method and in the analysis of the results we use an 
Ising-type spin model in order to elucidate the effects of mag- 
netic interactions on the energetics of the alloy. The rest of the 
paper is divided into two main sections and conclusions. The 
theoretical tool is briefly reviewed in Section II and the results 
are presented and discussed in Section III. 



II. METHOD 

The calculations are based on the density functional 
theory 36,37 and performed using the Exact Muffin-Tin Orbitals 
(EMTO) method3^ The EMTO method is an improved 
screened Korringa-Kohn-Rostoker method, 40 where the one- 
electron potential is represented by large overlapping muffin- 
tin potential spheres. By using overlapping spheres, one de- 
scribes more accurately the crystal potential, when compared 
to the conventional non overlapping muffin-tin approachi^i^ 3 - 

The EMTO basis set included s, p, d, and / orbitals. 
The one-electron equations were solved within the scalar- 
relativistic and soft-core approximation. The generalized gra- 
dient approximation in the PBE form was used for the ex- 
change correlation functional^! The EMTO Green's function 
was calculated self-consistently for 32 complex energy points 
distributed exponentially on a semi-circular contour, which 
included states within 1 Ry below the Fermi level. In the 



one-center expansion of the full charge density, we adopted 
an l-cutoff of 8 and the total energy was calculated using the 
full charge-density technique. 38,42 For each alloy the calcu- 
lated equilibrium lattice constant was used. The convergence 
of the total energy with respect to the number of k- vectors was 
tested. It was found that 1240 uniformly distributed k-vectors 
within the irreducible wedge of the Brillouin zone was enough 
for the present purposes. 

The Fe-rich Fe-Cr alloys adopt the body centered cubic 
(bcc) phase of a-Fe. A number of previous works demonstrate 
that below the magnetic transition temperature ( 900-1050 K) 
the energetics of Fe-Cr alloys with less than 10% Cr is well 
described using the substitutionally disordered ferromagnetic 
bcc phase. 4 ^ 6 - 

Since in the present investigation we map a large concen- 
tration region with Cr content approaching zero, the conven- 
tional supercell method would require enormously large su- 
percells. Here, we resolve this difficulty by employing the 
Coherent Potential Approximation (CPA). 47 Within the CPA, 
the alloy components are embedded in an effective medium, 
which is constructed in such a way that it represents, on the 
average, the scattering properties of the alloy. The EMTO 
approach in combination with the CPA has been applied suc- 
cessfully in the theoretical study of various structural and elec- 
tronic properties of alloys and compounds 38 demonstrating 
the accuracy and efficiency needed for the present investiga- 
tion. 

The EMTO results are analyzed by using a spin model! 8 , 
where the magnetic moments obtained from the EMTO calcu- 
lations are the input data and the output energy is used to cal- 
culate the effective chemical potential. This procedure allows 
a more transparent way to understand the effect of magnetism 
on the chemistry of the Fe-Cr-X alloys. 



III. RESULTS AND DISCUSSION 

A. Corrosion rate versus balance between bulk and surface 
chemical potentials 

Our previous investigations demonstrate that the chemical 
composition of the close packed surfaces of Fe-Cr alloys fol- 
lows closely the characteristic threshold-like behavior of the 
corrosion rate of ferritic stainless steelsi 2 ^ In dilute, corrodi- 
ble, Fei_ x Cr x alloys (x<0.05) the surfaces are exclusively 
covered by Fe, whereas the Cr-containing surfaces become 
favorable when the bulk Cr concentration exceeds 10 at.%, 
which is close to the generally accepted Cr threshold for cor- 
rosion resistant alloys. The discovered threshold point in the 
concentration of bulk Cr was found to be a consequence of 
the reversal of the relative magnitudes of the bulk and surface 
effective chemical potentials (/ipc — MCr)> which in turn re- 
flects the peculiar electronic and magnetic structure of Fe-rich 
Fe-Cr alloys, 45 ! 46 ! 48 ! 50 ^ 

In the present work we extend the previous investigations 
to questions to what extent the balance of the chemical poten- 
tials of Fe and Cr in bulk alloy can be tuned by doping Fe-Cr 
with typical minor alloying elements commonly used in com- 
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mercial steels. As a reference for the alloying effect, we use 
the aforementioned threshold value of bulk Cr which corre- 
sponds to the breaking of the pure Fe surfaces. We calculate 
the effective chemical potential in bulk ((/jfc — MCr) bulk ) for 
Fe-Cr-X alloys containing 5 at.% X (X=A1, Ti, V, Mn, Co, Ni, 
Mo) and compare the calculated (/ip c — A*Cr) bulk (Fig.Q} to 
the corresponding effective chemical potential at the surface 
of binary Fe-Cr ((/^Fc — /^Cr) surfacc X determined in our earlier 
work^iS 

In alloys for which the surface energy of the third compo- 
nent X is higher than that of Fe, the surfaces are expected to 
be covered mainly by Fe atoms. Therefore, in these cases the 
effective chemical potential /ipc — Mcr at the surface can be 
approximated by that of the binary Fe-Cr alloy and in these 
alloys the threshold value of bulk Cr corresponds to the on- 
set of the breaking of the pure Fe surface by Cr atoms. On 
the other hand, for those alloys where the surface energy of 
the third component X is lower than that of Fe (Al and Ni, 
Ref. l60h . the surfaces of the Fe-Cr-X alloys are expected to 
be covered mainly by X atoms. In these cases the lowering of 
the threshold value of bulk Cr corresponds to the increasing 
of the driving force for moving an Fe atom, compared to that 
of Cr, from the surfaces to the bulk of the alloy. 

By comparing the effective chemical potential /iFc — Mo in 
bulk and at the surface we are able to estimate, for Fe-Cr-X al- 
loys, the threshold concentrations of bulk Cr, which again are 
expected to be related to the onset of the corrosion resistance 
in Fe-Cr based alloys. The first-principles EMTO results for 
the threshold values (c t hr) of bulk Cr in Fei-^-o.osCr^Xo.os 
alloys are listed in Table [I] 



X 


Al 


Ti V 


Mn Co 


Ni Mo Cr 


cthr (at. 


7c) 


5 


4 1 


9 6 


3 1 8 



TABLE I: Theoretical results for the threshold values of bulk Cr cor- 
responding to the reversal of the directions of the Cr and Fe driving 
forces towards the surface in Fei-z-o.osCr^Xo.os alloys. The data 
is obtained by determining the intersection point of the calculated 
(MFc - /io) bulk (Fig.[D and (^ Fc - /i Cr ) surfacc . 2 The last column 
(Cr) shows the threshold value for Fe-Cr alloys, which agrees well 
with the experimentally observed onset point of corrosion protection 
in Fe-Cr. 



B. Chromium threshold as an indication of the energetics of 
the Fe-Cr based alloys 

Plotting the data shown in Table U as a function of the lat- 
tice parameter (a) of the alloy (Fig. [2j reveals an approximate 
decreasing trend of the Cr threshold value with the volume of 
the alloy. The observed general trend shown in the relation 
between the Cr threshold value and the volume of the alloy 
can be understood by considering the total energy of the alloy. 
The approximate connection between the threshold value and 
the volume can be elucidated as follows: 



(i) The increased volume is related to the decreased bond- 
ing between the atoms, and consequently to the in- 
creased total energy of the alloy. 

(ii) The increasing total energy of dilute Fe-Cr alloys de- 
creases the magnitude of the negative slope of the mix- 
ing enthalpy of Fei^Cr^ at x » 0. This is due 
to the characteristic shape of the mixing enthalpy of 
Fei-zCr-r. alloys at small x£ 

{Hi) The relation between the effective chemical potential of 
bulk Fe-Cr and the slope of the mixing enthalpy (AH) 
of the alloyS! 

3 AH 

(MFo - Mcr) bulk ~ — + constant (1) 

shows that the effective chemical potential at x 
consequently decreases. 

(iv) Finally, the decrease of (hf c — jUo) at small x shifts 
the threshold value c t hr to lower Cr contents. 

At this point it should be noted that vanadium deviates signifi- 
cantly from the general trend of c t i u with volume showing the 
largest reduction in the bulk Cr threshold with increasing vol- 
ume (Fig. |2). Therefore, the threshold value of the bulk Cr as 
a function of minor alloying element needs to be investigated 
in more detail. 



C. Establishing the role of magnetic interactions in the 
thermodynamics of Fe-Cr and its derivatives 

Because the phase stability of Fe-Cr alloys depends signif- 
icantly on the magnetic interactions between the constituent 
atoms, 2dS r 50 i 62 i 63 it is instructive to concentrate on the mag- 
netic properties in the present case too. In order to trace the 
role of the additional alloying elements (X) on the magnetic 
contribution to the bulk threshold of Cr in Fe-Cr-X alloys we 
analyze our ab initio results by using a spin Hamiltonian 

Hi, = ($ + Sj)a i( Tj/2 + (1 - SiSfcitTjfi. (2) 

A similar Hamiltonian was employed by Ackland (Ref. l48h 
in his study of the Fe-Cr systems. The total energy of the N 
atom system (-Etot) per atom then reads 

1 N 1 14 1 N 

^tot/atom =jfY,2^ His= N^ Hi =< Hi >1 

i—1 j—1 2—1 

(3) 

where the j summation runs over the nearest and next nearest 
neighbors of the site j (14 sites in the bcc lattice 48 ). Applying 
equation (J3j to binary Fe-Cr alloy, we adopt Sp c = — 1, Sci = 
1, and (jpo and cro are the magnetic moments of Fe and Cr 
atoms, respectively. This yields ferromagnetic Fe-Fe coupling 
and antiferromagnetic coupling of Cr to its neighbors. Within 
the CPA formalism, we obtain for the Fei^Cr^ alloy 

< Hi> = x < H C r > +(1 -x) < H Fc > 

= 7{x 2 Oq t + 2x(l - SE)<TCrffFe 

- (l-x) 2 a 2 Fc ). (4) 
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FIG. 1: (Color online) Results of the EMTO calculations. Lattice parameter (a), magnetic moment per atom (m, in Bohr magnetons /^b), 
effective chemical potential in bulk (/ip c — fiCi) and the shift of the Cr threshold Ac t h r of Fei-^-o.osCr^Xo.os where X is the alloying 
component (Al, Ti, V, Mn, Co, Ni, Mo). Continuous black curve shows the results for Fe-Cr. Acthr as a function of X content is a linear 
approximation based on 5 at.% X data. 



Since we are interested in the relative energy balance of Cr in 
bulk, we extend the above spin model to Fe-Cr-X by merging 
the X atoms with the Fe atoms to form an effective Fe matrix 
into which Cr atoms are dissolved. For individual spins a we 
use the calculated magnetic moments (Fig. [I). The spin of the 
effective Fe matrix is taken to be the average of the magnetic 
moments of the Fe and X atoms weighted by their concentra- 
tions. As Fig. |3] shows, the alloying of Fe-Cr changes both 
the curvature and the position of the minimum of the total 
energy (mixing enthalpy), particularly around the point of 5 
at.% Cr which is an important region with respect to the bulk 
Cr threshold. Because the bulk effective chemical potential 
(/iFc — Mcr) bulk is a key quantity in determining the bulk Cr 
threshold and since the effective potential is approximately 
the derivative of the total energy (Eqn. (Q])) we calculate the 
derivatives of the spin model total energies from Fig. [3] The 
result is shown in Fig. |4] The derivatives of the total en- 
ergy of the spin model have approximately the same shape 
as the EMTO effective chemical potentials (fip e — /icr) bulk , 
shown in Fig. Q] Plotting in Fig. [4] the surface chemical po- 
tential (^Fo — /io) surfacc of Fe-Cr alloy, we can estimate the 



bulk Cr thresholds Ac t hr by determining the positions of the 
intersections of the (/ip c — /icr) bulk and (^tp e — /Lto) s ace 
curves. These estimated critical compositions are shown by 
green crosses in Fig. [2] The spin model based Cr threshold de- 
viates significantly from that of the EMTO calculations in the 
case of V and Mo. What comes to the qualitative agreement, it 
is of course expected that the spin model cannot compete with 
the ab initio calculations, but what is more important here is 
that the general trend predicted by the spin model agrees with 
the EMTO result. This suggests that the spin model is capa- 
ble of capturing the main features of the Fe-Cr-X alloy and 
to reveal the underlying magnetic interactions. The obtained 
shifts of the thresholds of bulk Cr ( Acthr) due to the alloying 
in the spin model agree surprisingly well with those obtained 
directly by the EMTO calculation (Fig. |5j. Presenting the 
results using relative quantities, 

Act hr = cthr (Fe - Cr - X) - c thr (Fe-Cr) 
c thr (Fe - Cr) c thr (Fe - Cr) 

shown in Fig. [6] it is realized that there is a significant correla- 
tion between the relative thresholds of the EMTO calculation 
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FIG. 2: (Color online) Threshold of the Cr concentration in 
Fei-^-o.osCr^Xo.os alloys (X = Al, Ti, V, Mn, Co, Ni, Mo; base 
refers to Fe-Cr) as a function of the lattice parameter. The brown 
circles show the EMTO results and the green crosses show the spin- 
model results. The straight line is a guide for an eye showing the 
general trend with the volume. 
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FIG. 3: (Color online) The total energy per atom obtained from the 
spin model. The symbols refer to the results of the spin model (Eqn. 
J4]() and curves are functions fitted to the spin model data. The inset 
shows the same results but from 1 to 20 at.% Cr emphasizing the 
almost linear behavior of the energy as obtained from the spin model 
within the 10-20 at.% Cr range. 



and the spin model estimate demonstrating definitely that the 
effective chemical potential of Fe and Cr in the investigated 
alloys is to a large extent determined by magnetic interac- 
tions. We close our discussion by analyzing the EMTO results 
shown in Fig. Q] The general trend seen in Fig. Q] is that all 
investigated alloying components lower the (/ipc — A*Cr) buIk 
for low-Cr alloys. When Cr content increases, this trend is 
diminished or reversed (Mn). When Cr content reaches 15- 
20 at.% the effect of the additional alloying is small. Man- 
ganese shows the most peculiar behavior: below 5 at.% Cr, 
Mn strongly decreases (/ipc — A*Cr) bulk whereas above 10 at.% 
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FIG. 4: (Color online) The derivative of the spin model total energy 
per atom (spin model, Fig. [3j multiplied by —1 and shifted to the 
common zero value at Cr at.% = 20. The short gray line segment 
shows the level of the EMTO surface effective chemical potential 
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, from Ref. 2) relative to the Fe-Cr curve. 
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FIG. 5: (Color online) The shifts of the bulk Cr threshold Ac t h r 
from that of the Fe-Cr case (EMTO result: 8 at.%, spin model result: 
7 at.%) due to 5 % alloying with Al, Ti, V, Mn, Co, Ni, and Mo, ob- 
tained by the EMTO method (brown bars) and the spin model (green 
bars). 



Cr it increases (/ipc — MCr) bulk - 

Maybe the most interesting data in Fig.[T]is related to the 
magnetic moments of Mn. It provides an illustrative evidence 
for the decisive role of the magnetic moments of the alloy- 
ing elements in the relative balance between the bulk chemi- 
cal potentials of Fe and Cr in Fe-Cr-X alloys. While Cr con- 
centration changes from to 5 at.% the Mn moment changes 
approximately from —1 to +1 Bohr magneton (/j,b) i-e- prac- 
tically through the whole range of the moments of the inves- 
tigated additional alloying elements (from V to Co). Corre- 
spondingly the (/ipo — Hcr) hulyi changes from the value of the 
Fe-Cr-V case to that of the Fe-Cr-Co case revealing the cm- 
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fective chemical potential (/ip e — /j,c, 
the magnetic moment of the doping atom plays a significant 
role in the energy balance of Fe and Cr atoms in Fe-Cr-X al- 
loys. The employed computational method provides an effi- 
cient tool to scan the effects of different alloying components 
at any concentration level. Since the stainless property of Fe- 
Cr based alloys is related to the presence of Cr near the alloy 
surface, our findings are useful in understanding and explain- 
ing the corrosion resistance of the Fe-rich Fe-Cr based alloys. 
Using the obtained results for the bulk it is possible to esti- 
mate the effect of a third component on the relative content of 
Fe and Cr at the surface. 



FIG. 6: (Color online) Relative shift of the threshold Cr concentra- 
tion Ac t hr/cthr(Fe - Cr) in Fei-^-o.osCr^Xo.os with (X = Al, Ti, 
V, Mn, Co, Ni, and Mo; base refers to Fe-Cr). The line shows the 
linear fit to the data f(x) = a + bx, a = -0.028 b = 1.29. 



cial role of the magnetic moment of the alloying element on 
the chemical potentials in Fe-Cr alloys. 



IV. CONCLUSIONS 

Using density functional theory, implemented in the 
EMTO formalism, and a spin model, we investigated the ef- 
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